Jet Quenching via Gravitational Radiation in Thermal AdS 



Edward Shuryak, Ho-Ung Yee, and Ismail Zahed 

Department of Physics and Astronomy, State University of New York, Stony Brook, NY 11794 

(Dated: October 5, 2011) 

We argue that classical bulk gravitational radiation effects in AdS/CFT, previously ignored be- 
cause of their subleading nature in the 1 /Ac-expansion, are magnified by powers of large Lorentz 
factors 7 for ultrarelativistic jets, thereby dominating other forms of jet energy loss in holography 
at finite temperature. We make use of the induced gravitational self-force in thermal AdSs to es- 
timate its effects. In a thermal medium, relativist ic jets may loose most of their energy through 
longitudinal drag caused by the energy accumulated in their nearby field as they zip through the 
strongly coupled plasma. 

PACS numbers: 



I. INTRODUCTION 

This paper stems from two different motivations. One 
is a (century old) quest to understand the radiation re- 
action force in classical electrodynamics and general rel- 
ativity, especially in higher space-time dimensions. The 
second (practical) motivation is the need to estimate the 
magnitude of such a classical force on a falling gravi- 
tational body in thermal AdSs space-time. This latter 
problem is related to the empirically important issue of 
jet quenching in strongly coupled quark-gluon plasma as 
unravelled by recent collider experiments at RHIC and 
LHC. 

Because of these two different motivations, we make 
the introduction a bit large to allow for a brief descriptive 
of the history of both subjects which is perhaps needed 
for some diverse readers. A more thorough analysis of 
the radiation reaction force in electrodynamics in higher 
even-space+time dimensions will be presented in a com- 
panion paper pLj. 



A. Radiation reaction force 

One general theoretical problem, with which theorists 
such as Abraham, Lorentz, Dirac and others struggled 
from about the beginning of 20-th century, deals with 
the understanding of the radiation reaction force. The 
concept of non-relativistic and relativistic classical radi- 
ation loss can be found in most textbooks on classical 
electrodynamics. For instance, the relativistic form of 
the Larmor formulae for the power radiated by an accel- 
erated charge reads 

P = -— f -x, (1.1) 

with x^{r) being the particle trajectory parameterized 
by its proper time r. The double dots reflect the proper 
accelaration. It is not difficult to write down a longitu- 
dinal force t hat drags the charged particle matching the 
energy loss (|1.1[). This procedure refers to the radiation 



or large distance derivation of the local foce, based es- 
sentially on the energy-momentum flux radiated on an 
asymptotic surface at infinity (large sphere). However, 
there are many reasons for why a local derivation of the 
self-force that only makes use of the local fields acting on 
the accelarated charge is warranted. 

There is a large body of research on this approach 
for both scalar, vector (electromagnetism) and tensor 
(gravity) radiation both in flat and curved space-time. 
Once the energy loss is accounted for, the radiation 
reaction force is fixed by the condition of balancing loss. 
Putting aside some known paradoxes related to this 
particular construction, we would like to focus on the 
following specific questions: 

(i) Is the radiation reaction force amenable to a 
local definition from a regulated self-field at the instan- 
taneous position of the charge? 

(ii) In electrodynamics the radiation reaction force is 
expressed via the acceleration, but in general relativity 
particles move along geodesies for which the covariant 
acceleration is zero. So what makes up the force? 

(iii) The Larmor formula involves the acceleration 
thanks to the dipole approximation in non-relativistic 
electrodynamics, but what about higher derivatives in 
higher dimensional settings and general relativity? 

While we do not address all of these difficult issues 
here, we aim at exploring the practical importance of 
some theoretical proposals already in the literature. If 
(i) is true, it would simplify calculations considerably, 
especially in the case of general relativity since defining 
distant surfaces in curved backgrounds is in general elu- 
sive. More specifically, can the radiation reaction force 
be the self-interaction force following from the charge's 
own field due to its past trajectory? 

The closest to formulating an answer to a local defini- 
tion of the self-force in four-dimensional gravity was the 
work of Mino, Sasaki and Tanaka and also Queen and 
Wald [2j |3] , where the gravitational radiation reaction 
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force was obtained in a closed form 
mx'' = Gm^x^x'' / dr 



(1.2) 



whereby the integral is carried over the past world line 
of the particle up to a regulated time r_. Here G~ is 
the retarded propagator for the Einsten equation with 
the particle's own stress tensor acting as a source. While 
more details on this equation will be given below, here we 
just note that the bracket in (1.2) is just the Christoffel's 



force from a metric perturbation as induced by the past 
history of the particle. Note that the result (1.2) was 



derived in harmonic gravitational gauge: in this work we 
will ignore all subtleties that may be associated with this 
particular gauge choice. 

Does this force or its simplified analogues really work? 
Clearly, for cyclotron radiation in flat four dimensional 
space-time, a force given by an analogous expression 
would vanish. The simplest way to see this is to note that 
the retarded propagator in flat four-dimensional space- 
time reads 

G-(x,xO = ^e(xP - x''')^(ix - x'f), 

which is totally localized on the light cone. As a result, 
there are simply no points on the cyclotron world line 
that intersect the past light cone. However this draw- 
back is short lived if were to notice: 1) In curved space- 
times the retarded propagator gets contributions from 
light rays reflecting on the curved background thereby 
causing intersections with the past light cone; 2) In even- 
space+time dimensions the retarded propagator has sup- 
port on the entire interior of past light cone even in flat- 
space, causing contributions from the past history of the 
particle. (1.2) is a valid local radiative self-force in grav- 
ity. 

We will use ( |1.2[ ) to address the issue of jet quenching 
in AdS/CFT correspondence. Our main observation is 



that for ultrarelativistic jets, (1.2) is to be used in the 



regime of very short proper times such that the Taylor 
expansion of all the factors in the numerator in proper 
time would provide the flrst flnite term after all inflnities 
are subtracted and renormalized. In our analysis, the 
magnitude of the proper time e is 



1/7' « 1, 



(1.3) 



as will be made clear in later sections. In the small e- 
expansion the radiative self-force in curved and thermal 
AdSs follows from the flnite term in (1.2) after using the 
leading singularity of the retarded propagator G~ with 
the appropriate Taylor coefficient. 



B. Jet quenching 

Jets are the most spectacular aspect of QCD physics 
in high energy setting. Their interaction with the 



quark-gluon plasma in heavy ion collision has long been 
deemed very interesting, starting from the early work 
by Bjorken [4 . Indeed the phenomenon known as jet 
quenching has been observed in heavy ion collisions at 
RHIC; see early summary in Refs.[5, 6 . The observed 
magnitude of jet quenching was signiflcantly stronger 
than predicted by the perturbative estimates. 

More spectacular jet quenching phenomenon has been 
reported in the flrst LHC heavy ion run in November 
2010: events with huge energy (^ 100 GeV) having been 
lost and dissipated into the quark gluon plasma. Many 
surviving jets retaining a fraction of the initial energy 
and travelling back-to-back to the trigger jet appear to 
fragment outside the medium in the standard way. There 
is no indication of collinear radiation in the jet cone or 
nearby. The energy loss appears to be due to some longi- 
tudinal force applied to the moving charge as opposed to 
the expected transverse kicks due to perturbative QCD 
(pQCD). 

These observations strongly suggest the notion that 
the quark-gluon plasma is strongly coupled and that it 
needs to be treated in the AdS/CFT setting in trying to 
understand the mechanism of energy loss. Some of the 
previous works on this line will be briefly reviewed below. 
What seems particularly important is the dependence of 
the jet quenching on the jet energy E. For discussion 
purposes, let us assume that the energy loss is of some 
power p of the energy: 



dEjdx ^ EP. 



(1.4) 



In pQCD the energy loss through radiation results in 
splitting of the leading parton into many partons, each 
carrying a fraction of the leading energy. The scaling na- 
ture of pQCD implies that jet quenching is proportional 
to the energy itself; —dE/dx ^ E with p = 1. This is 
radically different from e.g. standard synchrotron radia- 
tion in accelerators which has a power p = 4. As noted 
by Pomeranchuck in 1939 and more recently by Kharzeev 
[8 , the energy loss with p > 1 leads to the following 



k{x)dx. (1-5) 



(p-i)Ef^ {p-i)Er 



The interesting result is a finite final energy Ef even 
if the initial jet energy is very large Ei ^ 00. If so, 
then all companion jets will have nearly same energy, 
independent of the energy of the original (trigger) jet! 

As we will show in this paper, the self-force from clas- 
sical gravitational radiation of a falling ultrarelativistic 
object in the holographic bulk does in fact produce a 
power p > 1 of the Lorentz factor. Thus this interest- 
ing behavior with higher energy jets stopping at shorter 
distances seems to be also a feature of the AdS/CFT 
correspondence. 
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C. Jet quenching in AdS/CFT 

The idea that strong couphng physics is at work in 
jet quenching and could be addressed in the context of 
AdS/CFT was initiahy put forward by one of us (with 
Sin) [9 . Since then, the number of apphcations and dis- 
cussions along these lines is legion. For completeness, we 
briefly review these discussions of jets in the AdS/CFT 
context. 

In QCD the simplest process in which jets have been 
first seen is e+e~ hadrons. The newborn quark- 
antiquark pair move in random direction, and create two 
back-to-back jets. Its description started from string 
breaking models at low energies (e.g. the so called Lund 
model) and later developed into an elaborate pQCD de- 
scription of the partonic coherent cascade at high ener- 
gies. The outcome is a good description of the data. 

Therefore it is not surprising that the first papers ad- 
dressing some of the issue in AdS/CFT set-up have the 
same setting: two charged particles moving with veloc- 
ity V back-to-back from each other [TOl [11] . (The same 
problem with non back-to-back velocities was addressed 
earlier in the context of the pomeron problem ^ ). As 
the charged particles recede in thermal AdS, a light-like 
string between them falls gravitationally in the bulk. The 
gravitational hologram on the boundary CFT shows no 
sign of a falling string, but rather a point-like explosion. 
Indeed, the boundary energy-momentum tensor T^^ in 
the comoving frame yields a spatial (pressure) part which 
is not isotropic. The falling string has nothing to do with 
a hydrodynamical explosion as no black-hole is formed. 
The strong coupling regime does not tolerate jets, even 
if the two charges move by straight lines as in e+e~ 
hadrons. More details about this point-like explosion, 
termed as conformal collider physics, was discussed by 
Hofman and Maldacena [13 . 

Ultrarelativistic jets in holography are well approxi- 
mated by light-like rays travelling near the UV bound- 
ary of AdS. Their demise comes at the hands of holo- 
graphic bulk gravity which cause them to ultimately fall 
towards the horizon and be absorbed by the standing 
black hole. This scenario was put forward originally 
in [9] , and recently revisited and sharpened by many [14]- 
[T6] (and references therein). Although a number of au- 
thors have considered gluons or quarks as end-of-the- 
string objects ^4] [15] falling on a static thermal black 
hole in AdS, it turns out that to assess the jet penetra- 
tion length it is simpler to approximate the end point by 
a path of near- null geodesic as initially suggested in [9 . 

The current issue on jet physics in the context of the 
AdS / CFT set up has now shifted to the question of what 
exactly are the masses of these objects and how the jet 
process is initiated. One of the main results is that the 
shape of the geodesic describes "scale evolution" of a jet, 
and is related to the jet quenching strength versus the 
length X travelled since its production. The geodesic path 
consists of two parts: a near horizontal one suggestive of a 
very nonlinear rise of energy loss —dE/dx ~ x^, followed 



by a near-vertical plunge into the horizon, finished by 
a jet explosion as it splashes on the black hole horizon. 
Needless to say, both results need to be carefully tested 
against the phenomenologically accessible information at 
both RHIC and LHC. We will postpone this to the future. 



D. Cyclotron versus gravitational radiation 

The first step toward relating two very different mo- 
tivations mentioned in the earlier part of introduction 
has been done by one of us (with Khriplovich) nearly 40 
years ago [17], applying the same method to 4 problems: 
cyclotron electromagnetic/gravitational radiations in fiat 
or curved 3+1 dimensional spaces in the ultrarelativistic 
regime 7 1. The results for the radiation intensity are 



rflat 



rflat 

grav 



I 



G4m'^-f^/R\ (1.6) 



grav 



where the Schwartzschield metric was used for the curved 
space. An ultrarelativistic particle was set to rotate at 
the (unstable) circular orbit of radius R ~ (3/2)r/j,. Note 
the decrease of the powers of 7 as we go from cyclotron 
radiation in flat space to gravitational radiation around 
a Schwartzchield black hole. Also note that the powers 
of energy in the formulae consistently exceed one. In 
the last case of (1.6), Newton's constant G4 appears in 
combination with the particle energy E = jm. As a 
result, if the particle is small and the massless limit m 
0,7 00, mj =drag fixed is considered, the limit is non- 
singular and thus the answer is the only possible one 
which makes sense physically. 

For pedagogical reasons, we would like to discuss the 
problem, i.e. gravitational radiation and self-force in 
AdS^ , in two steps. In this subsection we will compare 
the settings of the cyclotron radiation with that in pure 
or empty AdSs. Then in later section, we will introduce 
thermal AdS^^ and study our problem in two different co- 
ordinate frames complementarily; the plasma rest frame 
and the jet comoving frame. 

The familiar setting for cyclotron radiation consists of 
a point charge rotating on a circle of fixed radius R with 
7 = £^/m 1. The ratio of the proper time r to the 
coordinate time t is fixed, and e.g. 



X = 



dx"" 
~d^ 



7- 



dx"" 
Itt' 



(1.7) 



The 4-velocity is tangent to the world line, while the ac- 
celaration (induced by appropriate magnetic field) is or- 
thogonal to it and is directed radially toward the circle's 
center. All odd derivatives of coordinates share same 
direction with x^ while all even ones with x^ . An ex- 
tra power of gamma appears for each derivative over the 
proper time. While the cyclotron radiation is a textbook 
example in 3+1 dimensional space-time, it is not so in 
other dimensions. The difference, in respect to the self- 
force issue, will be discussed in 2+1 and 4+1 space-times 
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in the companion paper [T]. We note that these space- 
time dimensions are the relevant ones for holographic 
models of QCD in 1+1 and 3+1 dimensions. 



E. Motion in Schwartzschield versus AdSs 
space-times 

In the gravity setting the problem discussed in Ref. [17] 
is that of an ultrarelativistic particle rotating in (unsta- 
ble) circular orbit of a 3+1 dimensional Schwartzschield 
metric in polar coordinates (t,r, 6>,^). In contrast, an 
empty AdSs space is defined by the metric 



1 



-dt^ + {dx^)^ + {dx^f + {dx^f + dz^] (li 



which is a near- horizon limit of the 9+1 dimensional 
black brane solution. The brane itself is extended in 3 
coordinates 00 ^ 00 ^ 00 , and the metric is independent of 
them. These three coordinates are the analogue of the 
polar coordinates 0, (j) of the Schwartzschield metric, with 
z playing the role of the radial coordinate r. Unlike the 
Schwartzschield black hole, the empty AdS^ has no hori- 
zon since it originates from the so called extremal black 
hole with maximal possible charge. 

In the next section we will derive an expression for 
the self-force containing multiple derivativ es of the local 
trajectory over the proper time based on (1.2). To un- 



derstand its origin, we recall the geodesic equation for 
free falling in general gravity 
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g""^ {gcb,d + gdb,a - gcd,b) x'^x^, (1.9) 



where the right-hand side is driven by the Christoffel 
symbols. Since the Schwartzschield metric depends only 
on r, and the empty AdS^ metric only on the acceler- 
ation is directed only along these coordinates. The large 
tangent components of the 4-velocity on the right-hand 
side are of order 7^ and so is the accelaration. 

For the Schwartzschield case, equatorial geodesies fol- 
low by setting x^ = = 0. Thus 

-•^ = -^(1-7) i-"^iur + 2riu^f), (1.10) 

using X = u. The radius of interest for the ultrarelativis- 
tic circular orbit is a textbook result 



3r^ 



1 



1 



(1.11) 



Repeated differentiations of Tc yields further derivatives 
over proper time r if needed. Analogous expressions for 
AdSs can readily be obtained using similar arguments 
after setting the path in the geodesic inside the (t,x^, z) 
plane with x'^ = x^ = 0. 

In order to find the geodesic path itself, it is however 
more convenient to use existing integrals of motion rather 
than solving the second order differential equations for 



the acceleration just mentioned. The Schwartzschield 
metric does not depend on t,6>,0 and thus the lower- 
indexed momenta (per mass) UQ^ug^Ufj) are conserved: 
those can be called the energy and two components of 
the angular momentum, respectively. The AdS^ met- 
ric is independent of the 4 coordinates of the (bound- 
ary Minkowski world) (t, x^, x^, x^). As a result, all 4- 
momentum components Ua = Pa^ci = 1,2,3,4 are con- 
served. 

The geodesic line follows from the on-shell condition 
for the squared 4-velocity 



goo{u^) 



-gr 



RHS, 



(1.12) 



where RHS = —1 for a massive particle and for a 
massless one, and expressing all components in terms of 
the conserved momenta. For the Schwartzschield case 
one should use = g^^j, u"^ = g'^'^L and for AdS^, 



u''=g''''Pa. a = 1,2,3,4. 



(1.13) 



As a result, the respective radial velocities r or i can be 
expressed in terms of the coordinates. The result is a first 
order equation that is readily integrated. For instance, 
the trajectory in AdSs is 

dx^ z 1 



dx^ X P 



which integrates to the path 



= r dz' 

-J Zi 



(1.15) 



There is a qualitative difference between falling in a 
Schwartzschield background and in an AdSs background. 
Indeed, while in the former a trapped circular orbit is 
possible, it is not in the latter. Particles fall indefinitely 
in AdSs irrespective of their initial conditions. A jet in 
AdSs assumes a particle being very energetic with 7 1 
along x^ with 4- momentum conservation. The falling tra- 
jectory is composed of two parts: 1) an almost straight 
initial path with the x-momentum much larger than the 
z-momentum with X\ ^ i; 2) a final diving or almost 
vertical plunge into the horizon. 

Falling in the Schwartzschield and in the AdSs involves 
transverse spatial acceleration which is normal to the ve- 
locities (at least initially in AdSs). In the Schwartzschield 
case the circular fall produces a radiation intensity of or- 
der 7^ as shown in (1.6). Is the fall in the empty AdSs 
also accompanied by a 7^ radiation? The answer is neg- 
ative and one can readily see that there is no dependence 
of the radiation on 7. Indeed, since the AdSs metric is 
Lorentz invariant, we can switch to a frame where Pi = 0. 
The trajectory is then confined to the (t, z) plane with no 
transverse acceleration. This is not the case for thermal 
AdSs as we discuss in later sections. 
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II. SELF-FORCE IN GENERAL RELATIVITY 

The local self- force in 3+1 gravity with zero cosmolog- 
ical constant was derived originally by Mino, Sasaki and 
Tanaka and also Queen and Wald [21 13] . As we noted in 
the introduction and now we repeat for completeness, 



dr' 



1 



hca'h' 



(2.1) 



• a ' h 



with G being the graviton retarded propagator, 



□G 



aha'h' 



(2.2) 



where x = x(r), x^ = x(r'), J5(x,x') = d^{x — x')^/^ 
and 'g is De Witt's bilocal for parallel displacement along 
the geodesic ^ . Although the original derivation of 
(2.1) was carried in 3+1 dimensional space with zero 



cosmological constant and matter, its physical interpre- 
tation is applicable in any dimensions: the right-hand 
side is simply a modification of Christoffel symbols due 
to the retarded metric perturbation of the particle trajec- 
tory. Therefore we assume it to hold in general, especially 
for thermal AdSs in 4+1 dimensions. 

For ultrarelativistic jets, the eikonal limit is appropri- 
ate 



(2.3) 



-b^{x^x'\ Inserting (2.3) in (2.1) 



with □G-(x,x') 
yields 



mi^ ^ 47r2 G5 / dr' (2.4) 

^— CO 

((V^Q- -x"x^VdG-)(x-xO^ -2i;-i:VVVdG-) . 

In the above, we have dropped terms of the type V^^ 
as they are subleading in small proper time e-expansion 
than V^G". 

The scalar retarded propagator in a curved background 
of 4+1 dimensions can be related to the one in 2+1 di- 
mensions. Explicitly, 



G-(x,xO 



.+_f^ e(.',.)fc)^ 



, (2.5) 



where the expression inside the bracket is the retarded 
propagator in 2+1 dimensions. 0(x', x) is the generalized 
heaviside step- fun ction with a space-like surface through 
X (the final form (2.5) doesn't depend on the choice of 
this surface), and A is the Van-Vleck determinant 

A(x, x') = {g{x)g{x') f^ det (VaVa'<T(x, x')) , (2.6) 



which is a scalar two-point function, a is an another 
two-point scalar function 

a{x, x') = i(T - r') £ dr" x{t") ■ x{t"), (2.7) 

which is defined by the geodesic between x and x' . It is 
negative for time-like geodesies 



(j(x, x^) 



1 



(r-rO 



A2 



1 



d{x^ x' 



A2 



(2.8) 



where d{x^ x') is the chordal distance. While the latter 
is only defined locally for general curved space-times, it 
can be defined globally for dS and AdS spaces because of 
their spherical and hyperbolic nature. Indeed, for AdSs 
the finite distance is 



cos {d{x^ x')) — 1 



(x- 



2zz' 



(2.9) 



with z being the conformal direction and = —t^ + x'^ . 
This relation is readily derived by embedding AdSs in 
with a hyperbolic constraint. We further note that 
the r etarded propagator in AdSs is the known function of 
(2.5 ) derived in [19 , and our generic small time expansion 
is consistent with it. 



III. SELF-FORCE ON JETS IN ADS5 

A. The expansion in e 

For ultrarelativistic jets, the trajectory is character- 
ized by s mall prop er ti mes, and one can expand in it. 
Inserting (2.5) into (2.4) and following our arguments for 



the self- force in 4+1 dimensions show that the gravita- 
tional self-force for ultrarelativistic jets is dominated by 
the leading singularity in the (covariant) gradient of the 
propagator. 



V^G" 



3 VAct^ 
"4^ ' 



(3.1) 



with small e = (r — r') <C 1. The smallness of e will be 
explained further below. The problem is then reduced 
to a covariant expansion of for x' near x. For that 
we follow [iHl Ho] and expand a/A first covariantly in 
terms of the world function a and then proceed to Taylor 
expand a. Specifically 

VA = 1 + ^ R,5 a - ^ R,5;c ^ VV^ (3.2) 

288 360 60 

X(T a a a + • • • , 

For both empty and thermal AdS spaces there are sig- 
nificant simplifications, which come from the fact that 
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the Ricci tensor Hab is proportional to the metric, since 
the Einstein equation with cosmological constant tehs us 



1 _ A 



(3.3) 



with the scalar curvature R = —5 A/3. As a result, the 
covariant derivatives of the Ricci tensor simply vanish 



ab;c 



R, 



ah]cd 



0, 



which yields 



18 

1 

360 



648 

Rra n -Ti ^a^b^c^a 
a b^mcndO^ a a a 



(3.4) 



(3.5) 



after using the identity a^aa = 2a. 

For a falling jet in AdSs a tremendous simplification 
takes place when rewriting (3.5) in terms of jet trajec- 



tory proper time difference. Indeed, since the falling jet 
follows a geodesic at zero'th order, the bilocal a{x,x^) 
uniquely determines the distance crossed by the jet for 
short or small e. The geodesic definition and length 
are unique when is in the neighborhood of x. Thus 
a = and = ex^ where x^ is tangent to the jet 
geodesic. 

Another way to say this is to note that the covariant 
Taylor expansion of a along the particle trajectory (and 
not the graviton trajectory) to fourth order is 



12 



Dx • Dx, 



(3.6) 



with Dx^ = x^ -\- r^^x^x^, the covariant or long deriva- 
tive. The covariant derivative and its higher order vari- 
ants vanish along the falling jet geodesic making the first 
contribution exact. Thus along geodesic jet trajectories 
in AdSs 



360 



R" 



^Rr? 



ridjX X X X 



(3.7) 



neglecting renormalizations taking care of divergences 
an d sub leading terms in the large 7 -exp ansion. In terms 
of (3.7) the gravitational self-force (|2.4[) simplifies 



mx 



307r 



(/• 



de R" 



b^^mcnd X X X X X 



(3.8) 



The final integration along the proper time / de will 
be explained shortly, and is of order 1/7^. The drag- 
ging gravitational covariant self- force on the falling jet is 
longidudinal and of order 7^. It is entirely driven by a 
(squared) Riemann tensor. 



B. Final integration and magnitude of e 

So far we have assumed e, the proper time between the 
emission of the graviton and its action on the charge, to 



be small. We have not yet specified why and how small 
it is. A qualitative answer follows from the fact that 
the more relativistic the jet, the more localized its own 
trailing field. Indeed, the localization in time is expected 
to be t 1/7 and thus a localization in propertime r ^ 

A way to see this is to note that the final e integration 
needs to be cutoff at large e. Indeed, (3.6) in fiat space- 
time reduces to 



I flat 



12 



X ■ X 



(3.9) 



For ultrarelativistic motion, we can make the substitu- 
tion 



1/2 



e^e lH — e x ■ x 
V 6 

leading to a finite resummed result 



(3.10) 



de 1 



6 



-5/2 



V6 



X • X 



(3.11) 



with most of the contribution stemming from the range 
e ^ 1/7^. As a result, the dropped terms in our quasi- 
local analysis of the gravitational self-force are all sub- 
leading in 1/7. This argument shows how a schematic 
resummation of the subleading corrections yields a finite 
result for J de. The qualitative character of this subtitu- 
tion does not fix the overall coefficient exactly. For that, 
more quantitative work is needed. 

The scaling of e ~ 1/7^ is also in effect in the syn- 
chrotron analysis in even-space+time dimensions as we 
discuss in the companion paper [T] . This underlines again 
the analogy between cyclotron and gravitational radia- 
tion as we have noted above. 



C. Holographic gravitational force estimate 

The previous analysis is purely 5-dimensional, and 
to apply that to AdS/CFT, we have to translate 5- 
dimensional results into 4-dimensional gauge theory ones. 
The holographic duality holds in the double limit in 
which the number of colors is large compared to 't Hooft 
coupling Nc ^ X = Nc ^ 1 taken also to be large. In 
AdS/CFT 



27Vr 



(3.12) 



The string length Ig and the AdS radius L are tied by 
L'^/lg = A in the double holographic limit. For a mean- 
ingful gauge theory result, L which involves Ig should 
drop at the end, and we will show this explicitly in our 
result shortly. 



7 




-4.5- 



FIG. 1: (color online) Example of a geodesic trajectory for a 
falling particle in the (inverted) 5-th coordinate — z as a func- 
tion of x^. The open (black) circles are for massive particle, 
(blue) asteriks are for massless geodesic. The parameters are 
explained in the text 



The 5d metric has the form 



(3.13) 



where the gauge theory spacetime is measured by (t,x). 
From this, one can obtain relations between the local 5d 
energy EAdS and the 4d energy E as follows 



EAdS 



-E. 



(3.14) 



Similarly, the 5d (proper) time and the 4d (proper) time 
are related as 



L 

tAdS = —t. 

z 



(3.15) 



As a result, the 4- velocity x in the above result which is 
defined in terms of the 5d proper time is roughly 



L 



(3.16) 



where 7 is the ordinary Lorenz factor in 4-dimensions. 
Also, since m is the 5d rest energy it should be rescaled by 
z/L to tie with the 4d rest energy 777. 4. W e now note that 
the square of the Riemann tensor in (3.8) is of dimension 
but independent of L, 



R 



2 ^ P{z/zh) 



(3.17) 



with P a polynomial of z/zh ~ T/Qs. Qs is the jet satu- 
ration scale to be fixed below. The final integration J de 



is of order times some function of T /Qs. The 1/7^ 

was explained before, whereas the factor L follows the 
factor in ds^ ^ since the geodesic motion is indepen- 
dent of the overall factor of . 

With these considerations, the 5-dimensional non- 
covariant force which is equal to 5-dimensional gravita- 
tional energy-loss takes the form 



dEAdS 

dtAdS 



L3 



2 dE 

~dt 



de I {x) 



N2 {^4) VW (^^) ^^"'"'^ 

Y{z/zh). (3.18) 



/^\2 (m47)^ 



With this, one can make a safe massless limit by replacing 
= £J as it makes sense physically. Thus 



dE 
~dt 



E^ 

m 



(3.19) 



with F(0) =0. L drops in the final result for the bound- 
ary observables. We set L = 1 in the rest of our paper. 
In thermal AdSs the universal function F appears solely 
dependent on the ratio z/z^ with the holographic direc- 
tion playing the role of a renormalization scale evolution. 
Indeed, along the longitudinal direction, the jet forms 
at Zi ^ l/Qs and ends at Zh = l/(7rT) the black- hole 

horizon. 

Although ( |3.19 ) is subleading in l/A/'c, it is magnified 
by a large factor 7^ 1 thereby dwarfing this suppres- 
sion. Since the gravitational radiation effect is classical, 
this longitudinal dragging force is likely present at the 



boundary. We note that the work following from (3.19) 
yields a power or intensity 



E^ 



dz 



(3.20) 



in agreement with (1.6). We cannot completly trace the 



origin of this agreement. 

Finally, we note that the longitudinal covariant force 
following from a dragging colored string is of the order of 
7a/AT^ [14, 15 . The ratio of the longitudinal drag from 
gravitational radiation (or selfforce) to color is 



gravity — radiation — drag 
color — drag 



7 



10 



2. .4 



mVx 10 X 5 



,(3.21) 



not small, for typical jets at RHIC and LHC with 7 = 
10 — 100. Of course, the derivation given is perturbative, 
without back reaction explicitly included. (It means it is 
only formally valid for A/'c, A exceeding the realistic values 
ofNc = 3 and A ^ 25.) 

In summary, we argue that the gravitational radiation 
drag is comparable to the color drag at RHIC, and is 
perhaps the dominant mechanism for jet loss at the LHC. 
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(Other effects subleading in l/A^c? such as the fur- 
ther terms supergravity Lagrangian beyond the Einstein- 
Hilbert Lagrangian, may also in principle lead to higher 
power of 7: we have not yet inverstigated their effect.) 

IV. MATTER VERSUS WIND FRAME FOR 
THERMAL ADS5 

We now proceed to analyze in detail the radiative grav- 
itational force for thermal AdSs described by a black- 
brane with a horizon located at some value in the 
5-th coordinate. Holography asserts that this classical 
set up is dual to thermal A/'=4 supersymmeteric gauge 
theory at the boundary with temperature T = l/(7rz/i). 

A. The geodesies in the (plasma) matter rest frame 

To assess the effect of the thermal bath on the longitu- 
dinal induced gravitational self-force, we use in this sec- 
tion the standard thermal AdSs metric with plasma mat- 
ter at rest. In terms of the coordinates (t, x^, x^, x^, z) 
the metric is 

ds^ = ^ (^-fde + [dx'^f + {dx'f + [dx^f + , 

(4.1) 

with / = 1 — [z/ ZhY. In this frame the jet is represented 
by a point particle moving with some 4- velocity = u^. 
At the start, the jet moves ultra-relativistically in the x^ 
direction, so that u^^u^ = 7^ with 7 1, while is 
small. As it proceeds, its longitudinal velocity decreases 
while its velocity along the holographic direction z direc- 
tion increases, till finally it approaches the vicinity of the 
horizon. 

Using conservation of the 4-momentum of the particle 
one can readily derive the following expressions for the 
components of the 4-velocity 

i = z^E/f, x = z'^P, (4.2) 
z' = z\E'-fP')-z'f, (4.3) 

from which the path 

fz p 

= dz' ■ (4.4) 

Jz^ ^E^-fP^-f/z'^ 

In pure AdSs, / = 1 leading to the Lorentz invariant 
combination — in the integrand. This is no longer 
the case in a thermal medium with / ^ 1. 

An example of the corresponding trajectory in the x-z 
plane is shown in Fig{l]by the open circles. For compar- 
ison we also show a geodesic for massless particles (blue 
asterisks). Note that a jet here as a particle starts from 
the initial value of the holographic coordinate 2:^, taken 
to be 1, and ends dit z = Zh- In this example = 5. 



After selecting the jet gamma factor E = 40 and fixing 
the initial = 1, the value of the jet momentum follows 
from the on-shell condition as usual. What we see in 
Fig{l] is qualitatively similar to the standard trajectory 
of a stone thrown from a cliff. The trajectory turns down 
as the increasing z-momentum takes over the conserved 
X- momentum. 

The choice of Zi is a choice of the particle mass or jet 
scale (in the renormalization group sense) which is also 
referred to as the saturation scale Qg- Phenomenology 
sets it to be 1-2 GeV at RHIC and 2-3 GeV at LHC. 
With the units thus fixed, the jet energy in this exam- 
ple should correspond to 40-80 GeV at RHIC and 80-120 
GeV at LHC. The corresponding matter temperature is 
in this example T = Qs/Stt, which is lower than the tem- 
perature occuring in real hadronic matter. We selected 
it, however, in order the stopping distance (about 24/ Qg) 
be (for Qs = 1 GeV ) about 5 fm, or comparable to the 
size of heavy nuclei used in the experiments. 

The elapsed time (not shown) follows closely the dis- 
tance except for a divergence as z ^ z^. Any falling ob- 
ject freezes in the distant observer frame as it approaches 
the horizon. In jet experiments, the time is usually lim- 
ited by the matter freezeout time, thereby the falling jet 
remains close to z^ but never crosses it. 

The metric and thus all tensors are independent of 
gamma, and thus the counting of the power of gamma 
is straightforward. It is reduced to counting the powers 
of the 4-velocities. For example 

x^^^^,xV^7^ (4.5) 

so that X • X ^ 7^. Of course, one should check that 
one gets nonzero result for an ultrarelativistic particle. 
An important exception involves contractions with the 
the Ricci tensor which has a special property Rab ^ gab- 
Therefore its contraction HabX^x^ ^ 7^ and not 7^ as 
naively expected. 

Unlike the Ricci terms, the one with the squared Rie- 
mann tensor does not have a simplified structure and 
thus retains its naive order of magnitude ~ 7^. More 
specifically, we found 

KKUUUU = K^^^s g'^'^g'^^ ^f,au(3 u'^u^u^u^ (4.6) 

^4(1-^4/^4)2' 

where P is some lengthy polynomial of degree 4 in the 
velocities with coefficients depending on 7, z. Instead of 
presenting it as such, we show in Fig j2]its dependence on 
the longitudinal distance x^ travelled rather than z. As 
we see from the plot, this contribution to the self-force 
not only depend strongly on the jet energy 7 but varies 
rapidly with the longitudinal distance crossed by the jet. 
It takes extremely large values near the jet final stopping 
point. Clearly, the back reaction from this force on the 
travelling jet needs to be taken into account. 



10"^ 
10-1 

10-2 
10-3 




FIG. 2: Dependence of the combination of two Riemann 
tensor and four velocities normalized to the leading power 
of the jet Lorentz factor 7, namely HHunuu/j"^, versus the 
longitudinal distance (in the direction) travelled by the jet, 
for the numerical example displayed in Fig.l. The units, as 
explained in the text, are fixed by the particle mass or Qs- 



B. The "hot wind" frame 

The dependence of all parameters in question on the 
jet energy (Lorentz factor 7) is also instructive to study 
in an alternative frame, in which the longitudinal jet mo- 
mentum Px can be set to zero. Such boosted frames were 
used before, e.g. in the study of moving charmonium [2T] . 

First, it is instructive to see how the cyclotron elec- 
tromagnetic radiation is recovered in this frame. In the 
original frame all powers of 7 follow from the velocity, 
while in the boosted frame the velocities are 7^. The 
enhancement in this frame is transferred to the magnetic 
field which is magnified by 7. 

The same happens with the boosted metric. In empty 
AdSs the combination —dt'^ -\-dx^ turns to —dt''^ -\-dx^ 
after boosting, which is clearly Lorentz scalar. In the 
boosted frame, the jet is at rest with all 7's gone. Empty 
AdSs is Lorentz invariant with free particle motion being 
uneventful. 

At finite temperature, manifest Lorentz invariance in 
the metric is lost with 



dt' 



ldx^')\ 



(4.7) 



and dx^ = after the boost. Thus the metric gets 
enhanced by 7^ in this case. This amounts to effectively 
boosting the energy density. In the rest frame of the jet 
the effective temperature is enhanced 



T' = 



(4.8) 



Note that at large 7 the boosted metric in (4.7) has non- 
zero gi2 . There is "wind" in the x-direction. This has 
obvious consequences on the particle motion. While it is 
still true that there is conserved energy and momentum 



-100 



-60 
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FIG. 3: (color online) The same example of a geodesic trajec- 
tory for a falling particle as in the previous figure, but in the 
boosted "hot wind" frame. The points are in the (inverted) 
5-th coordinate — 2; as a function of x^, which is now also 
negative as the "hot wind" blows in the negative direction. 



given by the lower components uq = E^ui = P, the 
upper components are mixed by the non-diaginal metric. 
As a result, even if P = 0, one gets a nonzero velocity in 
the X direction since 



(4.9) 



The upper component metric is obtained in the standard 
way with 



dt 
dr 

dx 



^^M±£V^, ,4,10, 



^6 



^4 • 



The geodesic path is again derived from the on-shell con- 
dition QahU^u^ = — 1 after solving for the z- velocity 



z\2 



-v6„,2 77.2 



E Zy^Z 



(4.11) 

^2 



This is readily integrated. As usual, the ratio {u^ /u^Y 
eliminates the proper time and yields the falling trajec- 
tory x{z). 

Fig|3] shows the falling jet geodesic in the hot wind 
frame for the parameters used in the matter rest frame. 
The motion is a free fall along the z-direction, followed 
by a horizontal motion along the wind direction. Ini- 
tially, the particle is nearly standing and non-relativistic 
with ^ 1 and all other components zero. Then these 
components get larger with typically ^ 

In the initial time, all the 7 dependence is carried by 
the wind term in the metric ~ [zj Z}^^^^ . Note that for 
our two Figures, the numerical value of the wind term 
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is ~ 40^/5^ = 2.56 which is not very large compared to 
1. It is 10 times more for twice more energetic jets, we 
now detail how this combination is enetring the various 
contributions to the gravitational self- force in this frame. 

First, let's consider the particle acceleration along the 
geodesic 



-pa -b-c 



(4.12) 



Initially the velocities are of order 7^, with all powers of 
7 carried by the Christoffels. Typically, 



p5 _ 
J- 11 — 



(4.13) 



Even when the particle is at rest with = 1 (and all 
others zero) there is an acceleration of order 7^ along 
the z-direction. The calculation of the Riemann tensor 
contribution to the self-force is straightforward in this 
frame. To leading order, the result is 

11.^.69^" g'^^K^^aupu^u^u'^u^ « jHz^P', (4.14) 

where P' is another fourth order polynomial in the ve- 
locities. For 2;^, = 1 it is 

P'= +3z^{u'f + 3z^{u^f 

-12z%u*fu'' + 18z\u^ u*f 
-I2z^{u''fu*-Bz^{u''f 
-482^(u^ u^f +?,2z'^{u'fu^ 
-8z^(zi*)^ + 32z^(w*)V 
+42(u^ w*)2-28(w*)V 
+7(w*)^ - 28(u^)3u* + 7(w^)^ 

(4.15) 

It is non-zero even in the initial part of the trajectory 
when the particle is nearly at rest = 1, li^ = li^ = 0, 
as it has the {u^)^ term. However, as the wind accelerates 
the falling jet, the polynomial is seen to keep the same 
power of 7. The quantity plotted in Fig. 2 is of course a 
scalar, so it is same in both frames. 



V. SUMMARY AND DISCUSSION 

Starting with a companion paper we first mention 
its main result: the longitudinal selfforce in flat even- 
space+time 

2 + 1: {mx'')L^ ^ n^ ^^ 



4 + 1 : {mx'')i 



>/3 • X 

^2 



X • X 



lOV^ 



(5.1) 



for an electromagnetically charged particle undergoing 
circular motion in 2+1 and 4+1 dimensions respectively. 
Both of these induced forces match the expected radia- 
tion asymptotically. 



Our main result is that th e lead ing longitudinal gravi- 
tational selfforce (3.8) with (3.11), i.e. 



mx 



G^rr? ( 4 



307r V V6x • X 



^ "R , r^e-b-c-d -a 

e b^^TTicnd , 

(5.2) 



derived using similar arguments, starting from ( 1.2 ). The 
structure itself appears from the leading short distance 
singularity of the propagator in 4+1 dimensions. Its part 
in brackets, coming from resummation of the higher order 
terms, is schematic and we do not claim exactness of its 
coefficient, just the power of 7. The convergence of the 
expansion follows from e ~ 1/7^ which is expected in 
the ultrarelativistic limit. Indeed, at large 7 the fields 
are Lorentz contracted causing the m to be localized over 
a proper time of order 1/7^. In (3.11) we have shown 
how this cutoff in e is obtained through a qualitative 
resummation of the higher order terms. More work is 
required to fix this resummation quantitatively. 

It is interesting to note the fundamentally different na- 
ture of the approximations leading to the radiation/self- 
force in the non-relativistic versus the relativistic limit. 
In the former, the radiation proceeds through large wave- 
lengths or small frequencies via a multipole expansion at 
large distances away from the charge/mass. In the latter, 
the expansion is dominated by the short wavelengths or 
frequencies, that is by the field in the immediate vicinity 
of the charge/mass. In a way, the two approaches are ex- 
pected to be tied by the structures of energy-momentum 
conservation. 

To assess the magnitude of the gravitational self-force 
in AdSs we have used two different frames: the plasma 
matter rest frame and the jet rest frame. While scalar s 
and vectors are expected to transform according to the 
expected lore of relativity, the use of these two distinct 
frames shed light on the nature of the contributions to the 
self- force. The dominant contribution to the covariant 
longitudinal gravitational drag is of order 7^. 

The longitudinal self-force we have discussed is very 
different from those discussed so far in the context of 
holography. In particular, it has nothing to do with the 
color charge of the jet but rather with its energy. The 
massive amount of excitation left behind by the jet pro- 
duces backward gravitational drag or pull, which we have 
found to be growing dramatically with the jet energy. 
This effect is classical, large and physically well separated 
from the variety of other subleading effects one may think 
of. It is by no means subleading although in holography 
it is formally down by a power of Newton's gravitational 
constant G5 ^ 1/A^c • 

The present work is the first indication for the exis- 
tence of strong self-induced gravitational forces in ther- 
mal AdSs that have a power of the Lorentz factor 7 larger 
than one. Many of the terms with higher powers do ap- 
pear for arbitrary paths but are canceled out in the re- 
sult, mostly because on geodesies there is no covariant 
acceleration. The more involved problem of a dropping 
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string with a longitudinal color drag is likely to depart 
from the geodesic motion and therefore involves higher 
powers. This issue will be addressed elsewhere. 

It is amusing to note that the 7^ dependence of 
the non-covariant force and hence the intensity radi- 
ated by energy conservation is the same as in 3+1 
Schwartzschield case ^7\. We cannot trace the connec- 
tion except the fact that the combination (7771)^ = E'^ 
in the radiated intensity is the only one that is finite in 
the massless limit m ^ and 7 ^ 00 with E fixed. 

This paper is intended as a qualitative indication for 
the existence and importance of classical gravitational 
self-forces in studying jet quenching via AdS/CFT cor- 
respondence. We have not tried to assess their effects on 
jets trajectories in the bulk or their stopping distance on 
the boundary quantitatively. For realistic values of Nc 
and A, (3.8) appears to be somewhat too strong for ul- 
trarelativistic jets to be treated in perturbation theory. 
Therefore more work is needed to assertain the role of 
these forces in a reliable calculational framework. 

It is very challenging to try to understand the origin of 
this effect from the gauge theory side without recourse to 
holography and bulk gravity. As a parting comment, let 
us speculate that bulk gravitons and gravity are perhaps 
related to Pomerons and effective Pomeron field theory, 
which practitioners in the field have tried to construct 
for the past many decades. 



VI. APPENDIX 

The structure of the singular terms in even-space+time 
dimensions can be readily obtained using the same ar- 
guments as in odd-space+time dimensions. Indeed, the 
retarded propagator in arbitrary space-time dimensions 
follows from 

= -Q{x,x')lmGH{x,x')\a+io, (6.1) 

following DeWitt's prescription between the Feynman 
and Hadamard propagator. In general AdSs the singular 
part of the Hadamard propagator reads 



Gh{x,x') 



(-2a)3/2 



(6.2) 



where a is the world function bi-scalar (2.7). Both g{a) 
and w{a) are regular functions of a. g{a) is fixed by 
\I\Gh = in leading order 



2d^gaa^g{Ba-5) = 0, 
after using that a^aa = 2a. Since 

□cr = 5 - a'^da InA, 



(6.3) 
(6.4) 
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it follows that 

da\n{g/^/A) = 0, (6.5) 
which fixes g up to an overall normalization. 
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